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Abstract 

While measurements of the hyperfine structure of hydrogen-like atoms are traditionally regarded 
as test of bound-state QED, we assume that theoretical QED predictions are accurate and discuss 
the information about the electromagnetic structure of protons that could be extracted from the 
experimental values of the ground state hyperfine splitting in hydrogen and muonic hydrogen. 
Using recent theoretical results on the proton polarizability effects and the experimental hydrogen 
hyperfine splitting we obtain for the Zemach radius of the proton the value 1.040(16) fm. We 
compare it to the various theoretical estimates the uncertainty of which is shown to be larger that 
0.016 fm. This point of view gives quite convincing arguments in support of projects to measure 
the hyperfine splitting of muonic hydrogen. 
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I. INTRODUCTION 



The hyperfine splitting of the ground state of the hydrogen atom is among the most 
accurately measured quantities |l[ Q| : 

AEfJp^ = 1 420 405 751. 7667 ± 0.0009Hz, (1) 

The relative experimental uncertainty in Eq-Qdoes not exceed 10~^^. The theoretical predic- 
tions for AE^^^, based on QED, are less accurate. This is partly due to the computational 
difficulties which increase very fast for the higher order terms in the perturbative expansion 
in powers of a and (Za) 0, Q], and to the limited precision of the fundamental constants 
involved (the Rydberg known to 10~^^, electron to proton mass ratio and a known to parts 
of 10~^ [3]). The main uncertainties come, however, from the insufficient knowledge of the 
structure of protons. Because of this, the comparison of theoretical results with the exper- 
imental value of Eq. ^ can not test QED beyond the contribution of proton polarizability 
effects of the order of a few ppm. To perform a more precise test of QED, therefore, either 
additional information on the electromagnetic structure of protons should be used, or the 
comparison should be done between theoretical and experimental results on the hyperfine 
splitting in hydrogen-like bound states of charged point-like leptons. Muonium is most ap- 
propriate for the latter; indeed, the recent measurements of the hyperfine splitting of the 
ground state of muonium P| have been shown to agree with theory up to 0.5 10^^ (see ^ 
and references therein), so that the correctness of QED results about the hyperfine splitting 
of hydrogen-like atoms has been experimentally confirmed at least with the same preci- 
sion. As for the use of information on the electromagnetic structure of protons from other 
sources in attempts to reduce the theoretical uncertainty of AE^^^, until recently there 
was no theoretical progress in this direction, and the only realistic idea for years seemed 
to be to estimate the proton structure contributions from complementary measurements of 
the hyperfine splitting in muonic hydrogen ^ 3l- We are now going to critically analyze 
this idea in the context of recent theoretical results on the polarizability of protons 0, Q], 
of the development of new precision spectroscopy instrumentation Jlj^, and of our better 
understanding of the dominating proton structure contributions to the hyperfine splitting 
in hydrogen-like atoms. This will lead us to the alternative point of view to look at the 
hyperfine splitting measurements in hydrogen and muonic hydrogen as measurements of the 
Zemach radius of the proton by assuming that all QED predictions are credible. This way 



we shall obtain a "first experimental" value of the proton Zemach radius from the hyperfine 
splitting of hydrogen which may be compared to theoretical values based on different proton 
form factor fits. 



II. HYPERFINE SPLITTING OF HYDROGEN ATOM GROUND STATE 

To analyze the various sources of uncertainty in the theoretical value of AE'^^^ we put 



it in the traditional form 



ll|: 



AE^^^ = (l + + 6'"^) (2) 
where is the Fermi splitting 12| expressed in terms of the electron and proton masses 



me,mp and the dipole magnetic moment of the proton fip-. 

while and 6^^"^ are correction terms related to higher order QED effects and to proton 
electromagnetic structure due to strong interactions. (Eq. |21 is only correct in the leading 
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order, since higher orders QED and structure effects mix up 
0{a^) (without distinguishing a from {Za)), S^^^ is given by 

5QED ^ + 3 2 ^ ^2(i^g2 - ^) _ ^ loga(loga - log4 + ^) + 18.984 x - + . . . (4) 

where Oe is the anomalous magnetic moment of the electron. Note that the expression for 
^QED ^Qgg j^Q^ involve the mass ratio rrie/mp] all terms which depend on proton mass or 
come from strong interactions are included in 5^^^ . In turn, 5^^^ splits into a "static" part 
frigid accounts for the elastic electromagnetic form factors of the proton and can be 
calculated using data from elastic scattering experiments, a part 5^°^ that comes from the 
internal dynamics of the proton and could only be evaluated using data on inelastic processes 
with protons, and a part 5^^^ describing the strong interaction effects outside the proton, 
such as hadron vacuum polarization: 5^^"^ = _|_ ^poi _j_ ^hvp^ Two types of "static" proton 
structure corrections are incorporated in (5"sid^ associated with the spatial distribution of 
the charge and magnetic moment within the proton and with recoil effects, respectively: 
frigid _ ^Zemach _|_ ^recoii^ rjj^g former has bccu Calculated by Zemach and may be put in 
the form j^l: 

5" = / , "^'^2 ^2 (-GE{-p')GMi-p') -l]= -2am,pRp + 0{a') (5) 



(p2 + a'^mi y 



where ruep = rriemp/ {rrie + mp) and Rp is the first moment of the convolution of the proton 
charge and magnetic moment distributions, also known as Zemach radius of the proton. 
The recoil correction S^^'^°'^^ denotes the contribution of all terms which depend on the ratio 
me/rrip] for sake of sirnplicity we skip here the rather lengthy explicit expression of S'^'^'^°^^ 



nplici 

fl. 



which may be found in The proton polarizability correction 6^°^ and the hadron vacuum 
polarization correction have been evaluated recently using the available data on the 

proton polarized structure functions and on electron-positron annihilation into hadrons. 
Eq. 12 may now be put in a more detailed form: 

AE^^^ = E^ (1 + 6^^^ + 5^'='^^'^ + 5''"™" + 5P°^ + 6^^") (6) 

It is important for what follows to review the order of magnitude of the various terms in the 
right-hand side of Eq. Inland briefly discuss the uncertainty of each of them. The uncertainty 
of E^ is due to the uncertainty of rUe, rrip, fip and a and does not exceed 0.01 ppm 
is dominated by the anomalous magnetic moment of the electron; the uncertainty of the term 
comes from the uncalculated terms of relative order O(a^) and higher and is estimated not 
to exceed 0.001 ppm ^. (5'''sid jg dose to 40 ppm and references therein); its uncertainty 
- 2-3 ppm - is due to the uncertainty of the actual values of the proton radii and is of the 
order of 5%. The recoil term 5^"^™'' ^ 5.68 ppm ^ adds little uncertainty to S"^^'^. As for 
the proton polarizability correction 6^"^, until recently there existed only the upper limit 
^po\ ^ 4 ppm »4j, obtained by ascribing the whole discrepancy between AE^J^^ and AE^^ 
to the contribution from 6^°^. The present value of the proton polarizability correction, 
calculated on the ground of experimental data on the polarized structure function j^, is 
1.6 ± 0.6 ppm. The hadron vacuum polarization correction S^"'^ ~ 10~® is much too 
small. The overall uncertainty of AE^fjf^ is therefore of the order of 2-3 ppm and is entirely 
due to proton structure effects. All discussed quantities are summarized in the leftmost two 
columns of Table InJ the righmost two columns of the Table contain the numerical values of 
the corresponding quantities for muonic hydrogen and will be discussed in next section. 

As already pointed out, the challenging comparison of AE^^ and AE^^^ for hydrogen 
can not be regarded as a test of the QED contributions of order 0{a^) and higher 

because of the significant overall theoretical uncertainty. We adopt instead an alternative 
point of view to assume that the theoretical values of 5**^^°, S^^p and 5^°^ are accurate 
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TABLE I: Magnitude and uncertainty of the contributions to the hyperfine spHtting of the ground 
state of hydrogen and muonic hydrogen from various correction terms. 





Hyd 


rogen 


Muonic Hydrogen 




magnitude 


uncertainty 


magnitude uncertainty 




1420 MHz 


0.01 ppm 


182.443 meV 


0.1 ppm 




1.1610-3 


< 0.001 10-6 


1.1610-3 


10-6 


frigid 


39 10-6 


210-6 


7.510-3 


0.110-3 


^recoil 


610-6 


10-8 


1,710-3 


10-6 


^Pol 


1.410-6 


0.610-6 


0.4610-3 


0.0810-3 


^hvp 


10-8 


10-9 


0.0210-3 


0.00210-3 



and use the experimental data to determine the Zemach radius of the proton Rp as: 

Rp = - (aE«^V^^ - 1 - - (J--^™" - 5P°i - /(2mep«) (7) 

The above assumption is justifiable since all four correction terms are objects of QED, the 
only difference of 5*^^^ and 5^°' from the former two being that their evaluation requires the 
use of additional phenomenological information beyond first principles. From Eq. [7| we get 
the experimental value R„ = 1.040(16) fm, where the uncertainty ±0.016 fm comes from the 

The Zemach radius of the proton is defined in terms of an integral of the charge and 
magnetic form factors of the proton GE{k) and GE{k) over space-like transfer momenta 
k, k"^ = — (see Eq. or equivalently, by the first moment of the convolution of the charge 
and magnetic moment distributions P£;(r) and PM{f) in coordinate space jl^ : 

Rp = Jd'rrJ dh' pe{v - r') pM(r') (8) 

The directly observable quantity which is most sensitive to the Zemach radius of the proton 
is the hyperfine splitting of bound systems involving protons (compare to the Lamb shift 
which is related to the proton r.m.s. charge radius (iSl-) The experimental value of -Rp sets 
important restrictions on the theoretical models of proton electromagnetic structure and, in 
particular, on the parametrization of proton form factors, in terms of which the theoretical 
values are calculated. Consider as an example the values of Rp calculated from Eq. El using a 
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few popular approximations of the proton form factors. Numerical calculations give Rp=1.02 
fm for the dipole fit, and -Rp=1.067 fm for the fit of unfortunately, no information on 
the uncertainty of the parameters of the fit is available, and no conclusions could be made 
on the compatibility of these values with the "experimental" one. Both fits are consistent 
with /ipG'^(— k^)/G'A/(— k^) «i 1 for < 5 GeV^/c^. To account for the recent experimental 
results on the form factor ratio [2^; "we also evaluated Rp by using the Simon's fit for 
either the charge or magnetic form factor and expressing the other one using the relation 
fipGE{-k'^)/GMi-k^) = 1 - 0.13 (k^ - 0.04) |2o|, and got Rp=lMO fm and Rp=1.073 
fm respectively. Though preliminary, these estimates show that the current theoretical 
uncertainty of Rp significantly exceeds the experimental one, and that the experimental 
results on the proton Zemach radius may be used as a test for the quality of models of the 
proton in the limit of low transfer momenta. 



III. HYPERFINE SPLITTING OF MUONIC HYDROGEN ATOM GROUND 
STATE 



Muonic hydrogen is the only other hydrogen-like atom in which the hyperfine splitting 
of the ground state could be measured with high precision. Due to the large muon mass 
m^/irtf. 2 10^, the binding energy of the ground state of muonic hydrogen is of the order 
of 200 Ry, and the radius of the muon orbit is ~ ao/200 so that the energy levels of muonic 
hydrogen are more "sensitive" to the details of the proton structure than the levels of normal 
hydrogen. The expressions in Eqs. I2I3I and El for the hyperfine splitting of a hydrogen-like 
atom and the various contributions to it apply for muonic hydrogen as well; however, the 
energy scale and the relative size of the various terms differ significantly from the hydrogen 
case (see the rightmost two columns of Table ITTl) The Fermi splitting now is 183 meV that 
corresponds to a hyperfine transition wavelength of 6.1 /im. The explicit form of the higher 
order terms in ^Q^^ may differ from Eq. |3] since different momenta are expected to give the 
main contribution in loop integrations, while proton structure modifies the proton voerteces. 
We do not know of any published explicit expression of 5^^^ for the ground state of muonic 
hydrogen. To our opinion the lack of interest in the topic is due to the unclear perspectives 
of an experimental verification of the theoretical results. Indeed, the analogous coi 
to the Lamb shift of muonic hydrogen were evaluated with a very high accuracy 
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soon as the proposal for the experiment j23| was about to be put forward. We are therefore 
convinced that the uncertainty of 5*3™ may be brought down to 0.1 ppm if necessary. The 
evaluation of the recoil terms S'^^'^"^^ may not be that easy since the mass ratio m^/mp ~ 0.11 
is much larger compared to hydrogen, and terms of order 0{a{'m^/mp)"'), n = 2,3,4, are 
all expected to contribute by more than 1 ppm. Most appropriate in this case might be 
an essentially two-body approach based on the quasipotential equation [2^ l^l ; we assume 
that the theoretical uncertainty may be brought down below 10~^ this way, and put 1 ppm 
in the rightmost column of Table ITTl Since we do not know of any published result on (5'''^™'', 
we take as estimate of the magnitude of the recoil effects the formula for the leading recoil 
term in muonium of Ref. |^ and get 5''^'=™! ~ — {3a/ it) {m^/ nip) \n{m^/mp) ~ 1.7 10^^. The 
proton polarizability correction has been evaluated using the same methods as for hydrogen: 
^poi ^ (^4 g ± 0.8) X 10"^ [8]. The term 6^"^ describing the hadron vacuum polarization was 
shown to contribute by approximately 20 ppm [2^, unlike hydrogen where it does not give 
any considerable contribution in hydrogen. 

With the content of Table ITTl in mind, we are now ready to discuss the information that 
would be provided by measurements of the hyperfine splitting of the ground state of muonic 
hydrogen atoms. 

The proton structure correction (5'^*'' in muonic hydrogen is enhanced (compared to hy- 
drogen) by a factor of 2 10^. Therefore, a measurement of AE^^^ in (yU~p)is can not be a 
good test of QED since QED effects are overshadowed by the proton structure corrections. 
Further on, in both hydrogen and muonic hydrogen, the proton structure corrections 6^^^ is 
dominated by two independent terms: the Zemach term S'^^^^'^ and the polarizability term 
^poi_ While the Zemach term is directly related to a well defined physical parameter - the 
Zemach radius of the proton Rp (see Eq. El, 6^°^ is expressed in terms of the form factors 
and polarized structure functions of the proton in an indirect and case-dependent way and 
is not associated with a single parameter. Compared to hydrogen, both these terms scale 
approximately as (m^/mg). This all brings us to the conclusion that opposite to what was 
believed by some authors 0, 3], the measurements of AE^^^ in hydrogen and muonic hy- 
drogen atoms are not complementary in a sense which would let us extract the values of 
two universal parameters of the proton, characterizing its charge and magnetic distribution 
and polarizability. However, if assuming that all terms in the right-hand side of Eq. [7| are 
evaluated correctly by theory, these measurements may be regarded as repeated experimen- 
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tal determination of the Zemach radius of the proton. While the discutible point in this 
assumption is the credibihty ot the theoretical evaluation of 6^°^ (if neglecting 5"^^^), the 
repeated measurements of Rp in hydrogen and muonic hydrogen are the best way to verify 
it: compatible values of Rp extracted from the hyperfine splitting in hydrogen and muonic 
hydrogen will confirm the reliability of the theoretical values of 6^"^ and vice versa. 

The accuracy of -Rp depends on the uncertainty of 5^°'; a measurement of the hyperfine 
splitting of the ground state of muonic hydrogen based on the results of ^ would give the 
value of Rp accurate to 1%. This would be more precise than the value obtained in the 
previous section because of the smaller relative uncertainty of the theoretical uncertainty 
of 5^°' (see Table inj, however, things may change with the more refined theoretical results 
to come in the future. As already mentioned, such an accuracy would fairly allow to filter 
the numerous theoretical estimates of Rp and detect a deviation of Ge/Gm from 1 by 
distinguishing the values of Rp obtained with and without account of the JLab experimental 
results [2^. It would be preferable for this purpose to have the value of Rp accurate to 0.5% 
or better, that requires in turn that the theoretical uncertainty of 5^°' be brought below 
3 10^^ and that the experimental error of AE^Jp^ not exceed 30 ppm. 

The muonic hydrogen Lamb shift experiment, currently in progress at PSI js^, may 
provide at a later stage as a by-product the hyperfine splitting of the 2S'-state of the muonic 



hydrogen atom with a relative accuracy of the order of 0.5 10~^. According to 21[, the 
Zemach correction S'^^^^'^ for the 2S'-state is again ~ 8 10^^. Unless the (yet uncalculated) 
polarization correction 5^°' happens to be anomalously large in this case, this measurement 
would therefore provide the value of Rp with an accuracy of 5-10% which is below the 
expected accuracy of the measurements in the IS state. The significant improvement of 
the accuracy of the proton rms charge radius expected from PSI experiment will not help 
imcreasing the accuracy of Rp either because Rp is independent of the proton rms charge 
radius. In order to determine the Zemach radius of the proton from the hyperfine splitting of 
muonic hydrogen atoms in the Is or 2s states with accuracy 1% or better, the experimental 
error of AE^^ should not exceed 50 ppm - a requirement that is not met by the PSI Lamb 
shift experiment 2^. An alternative experimental method, based on the response of the 
muon transfer rate from hydrogen to oxygen to the population of the {fi~p)is para-state 



was proposed recently j2^. The method takes advantage of the recent progress in the 
development of tunable lasers in the far infrared range around 6.1 /xm Q]. The efficiency of 
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the method has been demonstrated by means of Monte Carlo simulations. The experimental 
error limits have not been discussed; the main source of experimental uncertainty is expected 
to be the Doppler broadening of the transition lines. 



IV. HYPERFINE SPLITTING OF MUONIC HYDROGEN MOLECULAR ION 

In conclusion, we would like to briefly outline an alternative possibility for determining 
the proton Zemach radius from a measurement of the hyperfine splitting of the muonic 
molecular ion pfi~p. The hydrogen muonic molecular ions (pfJ^'p) are formed in the excited 
ortho-state with orbital momentum J = 1 in collisions of muonic hydrogen atoms with 
hydrogen molecules Q|: 

(fi-p) + H2^ [{Pf^p)jp ee]* (9) 

The formation rate Xp^p is proportional to the hydrogen target number density ip: Xp^p = 
ip/ipo X 2.210^ s~^ [30]. At high densities Xp^p exceeds the muon weak decay rate Aq = 
0.45 10^ s'^ so that the muons spend most of their life time bound in muonic molecular ions. 
The spin-orbit and spin-spin interactions of the protons and the muon split the ortho-state 
with J = 1 into 5 hyperfine state labeled with the quantun numbers (sF) of the total spin 

nn 

s = Spi + Sp2 + and the total angular momentum F = s -|- J |2QI, Qll]; the para-state J = 
has no hyperfine structure. The energy separation between the ortho-levels with s=l/2 
and s=3/2 is ~138 meV, while the separation within each of these groups is two orders of 
magnitude smaller. 

A measurement of the splitting between the hyperfine states with s=l/2 and s=3/2 
might be based on the strong dependence of the rate of spontaneous ortho-para transitions 
Xop = Xop{s) on s: Aop(l/2) = 7.210^s""'^, Aop(3/2) = O.SlO^s""*^. Since only the ortho-states 
with s=l/2 are initially populated at high densities, the observable ortho-para transition 
rate is close to Aop(l/2). The idea of the experimental method would be to use a tunable laser 
to stimulate transitions from the s=l/2 to the s=3/2 hyperfine states, for which Aop(3/2) 
is an order of magnitude smaller, which would result in a resonance drop of the observable 
rate of ortho-para transitions. Though such an experiment would require the development of 
tunable far IR narrow band lasers in the range of 8.8 /im, it has some significant advantages 
compared with measurement in gaseous hydrogen at room tempertures - suppressed Doppler 
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broadening, small muon stopping volume, high energy density of the laser beam, high rate 
of laser- stimulated ortho-para transitions - which deserve a careful consideration in future. 
As for the theory, the accuracy of the currently available results on the hyperfine structure 
of {pfi^p)j=i is limited to the leading order Breit and Zemach corrections and to the muon 



anomalous magnetic moment 



221 ; therefore, the theoretical uncertainty is currently 



about 10 ^. In order to determine Rp to 1% from pfi p hyperfine splitting measurements, 
next order QED effects have to be taken into account together with proton polarizability and 



molecular ion finite size effects 



mm 

so that the theoretical uncertainty be reduced by an 



order of magnitude to meet an experimental error below 50 ppm. This is a challenging task 
since {pfi~p) is a bound system of particles with comparable masses in which three-body 
relativistic dynamics may show up yet in the next-to-leading order. 



V. CONCLUSIONS 



By assuming that the numerical results for the various terms in the theoretical expression 
for the hyperfine splitting of the ground state of the hydrogen atom are correct (within the 
limits of the claimed accuracy), we have determined the value of the proton Zemach radius 
Rp. We have also demonstrated that comparison of the experimental value of Rp with theo- 
retical calculation is a sensitive test of the quality of the approximation of the proton form 
factors at low momentum transfer and of the possible deviation of the ratio of the charge 
and magnetic form factors from 1. To verify the credibility of the theoretical evaluation 
of proton polarizability effects, we consider the possibility of measuring the hyperfine split- 
ting of the ground state of muonic hydrogen atoms or, in a remote perspective, of muonic 
hydrogen molecular ions. 
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